Evolving properties of two dimensional materials, from graphene to graphite 



Mattias Klintcnberg,^'El Sebastien Lebegue,^ Carlos Ortiz, ^ Biplab Sanyal/ and OUe Eriksson^ 

^Department of Physics and Materials Science, Uppsala University, Box 530, 751 21 Uppsala, Sweden 
^ Laboratoire de Cristallographie, Resonance Magnetique et Modelisations (CRM2, 
UMR CNRS 7036) Institut Jean Barriol, Nancy Umversite BP 239, 
Boulevard des Aiguillettes 54506 Vandoeuvre-les-Nancy, France 
(Dated: March 13, 2009) 

We have studied theoretically, using density functional theory, several materials properties when 
going from one C layer in graphene to two and three graphene layers and on to graphite. The 
properties we have focused on are the elastic constants, electronic structure (energy bands and 
density of states), and the dielectric properties. For any of the properties we have investigated the 
modification due to an increase in the number of graphene layers is within a few percent. Our results 
are in agreement with the analysis presented recently by Kopelevich and Esquinazi (unpublished). 



INTRODUCTION 

The recent explosion of scientific activity around the 
newly discovered two-dimensional material, graphene, is 
unprecedented since the discovery of the high tempera- 
ture superconductors in the late 1980'ies. 
The uniqueness of this material, and the technological 
advantages it promises, gathers researchers from differ- 
ent scientific fields. As in many previous major scientific 
breakthroughs the main component for the success story 
of graphene was the actual synthesis of the material[6[. 
Unexpectedly Novoselov et al. were able to fabricate a 
truly two-dimensional material, and free standing single 
graphene layers, the building block of graphite, produced 
by means of exfoliation. Graphene can also be grown on 
SiC by epitaxial growth. In addition growth of graphene 
on catalytic surfaces (e.g. Ni or Pt) has been demon- 
strated. An insulating thicker material can be grown on 
top and after chemically removing the primary layer, one 
is left with a single atomic layer of graphene on an insu- 
lating substrate[l|. 

Research on graphene was initially motivated by a 
highly spectacular phenomenon, namely mass-less Dirac 
fermions. Although there is nothing relativistic in the 
single electron Kohn-Sham Hamiltonian describing the 
electronic structure of graphene, the band dispersion 
turns out to have a very unique property, at least 
from calculations based on the local density approxi- 
mation (LDA) and the generalized gradient approxima- 
tion (GGA) , which enables a comparison to quasiparticle 
states obtained from the Dirac equation, with an effective 
speed of light of 10^ m/s U, Q and zero rest mass. 

Studies of the electronic properties of graphene have 
revealed an ambipolar electric field effect, with high con- 
centrations and high mobility (up to 15.000 cm^/Vs). In 
addition, conductivity properties reveal ballistic trans- 
port on the submicrometer scale, which is unexpected 
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not least due to that the measurements are not done in 
ultra high vacuum, and hence many different molecular 
species are expected to be absorbed and act as scatter- 
ing centers. Furthermore, graphene is the only material 
know to date with a quantum Hall effect (QHE) at room 
temperature [3, and the QHE was shown to be somewhat 
anomalous in nature 0. The realization of a material 
with negative index of refraction for electrons has also 
been demonstrated in graphene, which is due to that the 
electron states in the valence band have a group velocity 
antiparallel to the k-vector. Hence, the Veselago lens, 
which has the unique property of having a resolution not 
determined by the wavelength, has been demonstrated in 
applications with graphene [sj. 

The focus of the present paper is to investigate how 
the electronic properties of graphene evolves to those of 
graphite, by a systematic theoretical study of 1, 2 and 3 
layers of graphene. We will make comparisons between 
the calculated results and existing data for graphite, and 
in this way we will shine light on how the electronic struc- 
ture of sp^ bonded C layers evolve from that of graphene 
to that of graphite. The properties we focus on here are 
the electronic structure, the dielectric function, e(cj), and 
the elastic constants. 



DETAILS OF CALCULATIONS 

The electronic structure and elastic constats have been 
calculated using a highly accurate fuUpotential linear 
mufhn-tin orbital (FP-LMTO) method ^ within the lo- 
cal density approximation (LDA). A 36 x 36 x 6 k-point 
grid gave converged results. The FP-LMTO and PAW 
(see below) results are found to agree nicely. A 2sp3sp ba- 
sis set was necessary to obtain a correct electronic struc- 
ture, c.f. the minimal 2sp basis set, and the muffin-tin 
radii were optimized to cover 90% of the nearest neigbour 
distance. 

The calculations of the optical properties have been 
performed by using the VASP (Vienna Ab-initio Sim- 
ulation Package) code [13, Hi, implementing the PAW 
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formalism [1^. We have used the PBE exchange- 
correlation functional[l^. To insure negligible interac- 
tion between periodic images, a large value (20 A) of 
the cell parameter 'c' was used. The convergence of the 
dielectric function is obtained by using a 80 x 80 x 5 
Monkhorst-Pack mesh IJ] . For the plane wave expansion 
of the wave function a 400 eV cut-off was used. 



RESULTS 
Electronic structure 

The energy bands for 1, 2 and 3 graphene layers are 
shown in Fig.l. For one graphene layer we find an elec- 
tronic structure which is similar to that found by others, 
for instance in Ref. 1^ 1^, 17 1. At the K-point two en- 
ergy bands of p^ character cut the Fermi level (E^), and 
the energy dispersion is (close to) linear with respect to 
the crystal momentum. This represents the bands ref- 
ered to as mass-less Dirac Fermion states. For two layers 
the number of energy bands doubles, and there are four 
sets of p^ derived bands close to the K-point. Due to 
the interaction between the graphene layers these bands 
split apart so that only two bands cut Ep, and the en- 
ergy dispersion deviates more from linear compared to 
the situation for one layer. For three layers a set of six 
Pz derived bands can be found close to at the K- 
point. None of the bands cut Ef, i.e. a small gap is 
introduced between four of these bands, whereas two are 
split, and are found further away either above or below 
Fip. In general Fig.l suggests that the more carbon lay- 
ers that are introduced the wider energy range does the 
set of p^ bands span. This saturates for bulk graphite 
where in addition to the degenerate bands at E^^ there 
is one band ~ 0.7 cV above Ep and one band ~ 0.7 eV 



below Ef.[18| 

The corresponding density of states are shown in Fig. 

m 

The three layers bandstructure look like a combination 
of the bandstructure of one graphene layer and two layers 
with the exeption that a small bad gap has been intro- 
duced (0.01 eV). We also note that the FP-LMTO resuhs 
presented here do not show the band-overlaps around the 
K-point at Ei? as do the tight-binding results in. [lit 

Some experiments, using either the de Haas van 
Alphen effect pH . or the quantum Hall effect [2^ . have 
shown that electrons and holes with linear dispersion 
relations could exist not only in graphene but also in 
graphite. In fact, it was shown using a tight-binding 
description of the electronic structure psj that Dirac 
fermions are existing in graphene multilayers if the num- 
ber of layers is odd (systems with a mirror inversion 
plane). However, tight-binding calculations depend on 
a given set of parameters, and although they can provide 
clues about the general mechanisms, they are not as pre- 
cise as ab-initio calculations. Therefore it was important 
to check if the above assumptions are still valid when 
obtained from ab-initio calculations. We computed pre- 
cisely the electronic structure of each system around the 
K point, as shown in Fig.l and found that indeed linear 
bands are present for one and three layers of graphene, 
but not for two and four. In particular, the linear bands 
in the three layers system show a slope very close to the 
one of graphene, suggesting that their Fermi velocities 
will be similar. 



Dielectric response 

The imaginary part of the optical dielectric function is 
calculated as (see e.g. RefllSl) 



el'{uj) oc — ^ < kn|pj|kn' >< kn'|pj|kn > x /k„(l - hn')S{ei^„' - ek„ - huj) (1) 

knn' 
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where < kn|pi|kn' > is the expectation value of the mo- 
mentum operator between band states \n > and |n' > for 
states with the crystal momentum k, i (and j)= x, y or 
z, and V is the volume of the unit cell of the crystal. In a 
calculation of graphene, which is a two dimensional ob- 
ject, one has to make a somewhat arbitrary choice of the 
volume of the C atoms of the two dimensional unit cell. 
We have chosen to use a volume of these C atoms which 
is the same as the volume of C in graphite. In practice 
then the calculations were made for a cell with extended 
c-axis, in order to simulate isolated C layers, and the 



numerical value of ^2 (w) was scaled to correspond to a 
volume V which is that of graphite. 

Our computed dielectric functions (real and imaginary 
parts, albeit without a Drude component) for graphite 
and for one, two, and three graphene layers are presented 
respectively in Fig. [3] Results for the x-component of the 
momentum operator (see Eq. 1) {exx) are presented in 
full (black) lines, while the results corresponding to the z- 
component of the momentum operator (czz) are in dashed 
(red) lines. The data in Fig. [3] agree rather well with 
published data for graphite [l3|. In Reffisl a comparison 
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Figure 1: The LDA electronic bandstructure for 1-, 2-, 3-layers graphene and graphite. The hnear bands around the K-point 
near the fermi level can be seen for 1-, and 3-layers, i.e. odd number of layers. 
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Figure 2: The density of states for 1-, 2-, 3-layers graphene and graphite. Colored lines show /-resolved partial DOS, s (dot, red), 
p (dash, blue), and d (dot-dash, green). Black full line shows Total DOS. 



4 



between experimental and theoretical data for graphite 
was made, and it was observed that the agreement was 
rather satisfactory. Also, our results a gree well with the 
one of MarinoDOulous|20l| . For graphite [l^l, the observed 
features are mostly due to transitions between tt and tt* 
states (for the 4 eV peak), and to transitions between a 
and a* states on the high-symmetry line between F and 
M. 

The data in Fig. [3]suggest that concerning the real and 
imaginary parts of e, the calculated values are rather in- 
dependent on the number of graphene layers, i.e. the 
curves in Fig. [3] are essentially independent on the num- 
ber of C layers. The main effect of the thickness can be 
found for , where for one single layer of graphene. the 
transition between tt and tt* states are forbidden [l8|, so 
tzz is exactly zero between and ~ 10 eV. For two and 
three graphene layers, these transitions are not stricly 
forbidden but they remain very weak. 

However, the imaginary part of the dielectric function 
of graphite is slighlty different for u> close to zero from 
the one of multilayers of graphene because the bands near 
the high-symmetry K ^ which are responsible for the tran- 
sitions at suck energies, are more splitted in the case of 
graphite. Therefore, the peak for w ^ of Imt^x is en- 
larged, and then ReCxx comes out differently for graphite. 



Elastic constants 

The theory of elasticity of three dimensional objects 
can be cast in a simple equation 



Table I: Elastic constants of graphite [in TPa] and of graphene 
[in Pa m] for 1, 2 and 3 C layers. The tension coefficients 
(711, 712 and 766) for graphene are given in 10* dyn/cm. For 



E{y, S) aE(v„,(i) + v„Y. 



i](:iSi(^j6j, 
(2) 

where E{Vo, 0) is the total energy of the undistorted sys- 
tem at volume Vb, the sums run over Voigt index 1-6, q 
takes the value 1 if the Voigt index is 1,2 or 3, and it takes 
the value 2 if the Voigt index takes the values 4,5 or 6. 
Furthermore, Ti is an element of the stress tensor and Cy 
is the elastic constant. 2J] For a two dimensional object 
like graphene the theory of elasticity becomes somewhat 
modified, as discussed e.g. by Behroozi [2^. Hence, the 
expression in Eqn.l is modified to, 

E{A, S) « EiAo, 0)+AoY, t,S, + C^i^^^^^^i^j^ 

i ij 

(3) 

where Aq is the area of the unit cell. For a three dimen- 
sional hexagonal lattice there are 5 elastic constants Cn, 
C12, Ci3, C33, and C55, which for the two dimensional 
hexagonal lattice of graphene reduces to only Cn and 
C12. Because the expression in Eqn.2 involves an area 
instead of a volume in front of the summation, the unit 
of the two-dimensional elastic constant is different than 
that of a three dimensional clastic constant, where the 



completness C66 - 


= (cii - ci2)/2 is 


also given. 




Material 


Cll (711) 


C12 (712) 


C66 (766) 


graphite (calc.) 


1.098 


0.154 


0.472 


graphite (calc") 


1.211 


0.276 


0.468 


graphite (exp.) 


1.060\ 1.109" 0.180^ 0.139'= 0.442^ 0.485" 


1-layer (calc.) 


358 (35.8) 


55.0 (5.50) 


152 (15.2) 


1-layer (calc") 


(40.6) 


(9.2) 


(15.7) 


2-layers (calc.) 


368 (36.8) 


47.3 (4.73) 


160 (16.0) 


3-layers (calc.) 


358 (35.8) 


54.5 (5.45) 


152 (15.2) 




unit is Pa. Hence the unit of the two dimensional elastic 
constant is m-Pa. We will below report on our calculated 
elastic constants of graphene in this unit. 

In Figs. [4] and [5] we show the calculated total energy 
versus distortions corresponding to the elastic constants 
Cll and C12, respectively. It should be noted that we 
display the energy per C atom. As is obvious from the 
two figures the different systems react very similarly to 
distortion with roughly the same energy cost. Hence, 
the expansion coefficients for bulk graphite as well as for 
graphene, bi-layer graphene and tri-layer graphene are 
all very similar. This is consistent with the fact that the 
chemical binding which is relevant for these two distor- 
tions is governed by the sp^ bonds, which are very similar 
for the four systems shown in Figs [3] and O . 

For graphite the calculated distortions correspond to 
a value of Cii= 1.098 MBar and Ci2=0.154 MBar. Both 
values reproduce with acceptable accuracy the experi- 
mental values, see tabid. This gives credit to the ac- 
curacy of the calculations and enables us to trust the 
elastic constants of graphene, bi-layer graphene and tri- 
layer graphene. The elastic constants for these systems 
are also listed in Table. I. Unfortunately we are not aware 
of experiental data with which to compare these num- 
bers, and hence our theory serves as a prediction. We 
note however that experimental studies of elasticity of 
graphene have been published recetly, reporting on the 
Youngs modulus. [2^ Michel and Verberck [i^] have cal- 
culated the elastic constants for graphite and tension 
coefficients (cy « 2"fij/c) for graphene using the Born 
long wave method to obtain phonon dispersion. Table. I. 
shows that the elastic constants for graphite reported in 
2?! are too large compared to experiment (10% and 70% 
for Cll and C12, respectively). Because we consistently 
report smaller elastic constant and tension coefficients 
compared to [l^ we are confident of the accuracy of the 
first-principles calculations in this paper. 
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Figure 3: Real (left column) and imaginary (right column) part of the dielectric function of graphite (upper panel), and of 
three (second upper panel), two (second lower panel), and one (lower panel) layers of graphene. 
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Figure 4: The calculated total energy as a function of distortion. Cn + C12 is calculated using a.^^ = a, (Xyy — a and with the 
other elements of the distortion matrix aij =0. 




Figure 5: The calculated total energy as a function of distortion. Cn — C12 is calculated using a^x = a, ctyy = —a and with 
the other elements of the distortion matrix = 0. 
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CONCLUSION 

111 this paper we have studied theoreticaUy several 
materials properties when going from one C layer in 
graphene to two and three graphene layers and on to 
graphite. The properties we have focused on are the elas- 
tic constants, electronic structure (energy bands and den- 
sity of states), and the dielectric properties. In general 
we find very similar behaviour for all studied systems. 
For any of the properties we have looked at the mod- 
ification due to an increase in the number of graphene 
layers is within a few percent. The largest effect due to 
the thickness is found for t^z which is zero in the energy 
interval of - ~ 10 eV for monolayer graphene, and non- 
zero for thicker layers, including graphite. The mililarity 
in elastic constants, Cn and C12, for the here studied 
systems is naturally due to that these constants are de- 
termined by the covalent in-plane sp^ hybrids, which are 
essentially the same and independent on thickness. Our 
results are in agreement with the analysis presented re- 
cently by Kopelevich and Esquinazi. [Slj 
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